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Abstract 

The  dot  diffusion  method  for  digital  halftoning  has 
the  advantage  of  parallelism  unlike  the  error  diffusion 
method.  The  method  was  recently  improved  by  op¬ 
timization  of  the  so-called  class  matrix  so  that  the 
resulting  halftones  are  comparable  to  the  error  dif¬ 
fused  halftones.  In  this  paper  we  will  first  review  the 
dot  diffusion  method.  Previously,  8^  class  matrices 
were  used  for  dot  diffusion  method.  A  problem  with 
this  size  of  class  matrix  is  that  enhancement  of  images 
is  necessary  before  halftoning.  However,  enhancement 
may  not  be  desirable  in  some  applications.  In  order  to 
eliminate  the  enhancement  step,  we  increase  the  size 
of  the  class  matrix  to  16^  and  optimize  the  class  ma¬ 
trix  for  a  set  of  gray  levels.  In  the  optimization,  the 
Human  Visual  System  is  used  in  the  cost  fimction. 
The  optimization  is  done  with  the  pairwise  exchange 
algorithm.  Since  we  increase  the  size  of  the  class  ma¬ 
trix,  we  are  compromising  the  parallelism,  i.e.,  the  al¬ 
gorithm  will  terminate  in  16^  steps  rather  than  8x8 
steps.  This  is  the  price  peud  for  avoiding  the  enhance¬ 
ment  step. 

Ke3nvords:  Halftoning,  dot  diffusion,  parallelism. 

1.  Introduction 

Digital  halftoning  is  the  rendition  of  continuous-tone 
pictures  on  displays  that  are  capable  of  producing  only 
two  levels.  There  are  many  good  methods  for  digital 
halftoning[5].  Ordered  dithering  is  a  thresholding  of 
the  continuous-tone  image  with  a  spatially  periodic 
screen.  In  error  diffusion,  the  error  is  ‘diffused’  to 
the  unprocessed  neighbor  points.  The  dot  diffusion 
method  for  halftoning  introduced  by  Knuth[4]  is  an 
attractive  method  which  attempts  to  retain  the  good 
features  of  error  diffusion  while  offering  substantial 
parallelism.  The  method  was  improved  by  optimiza¬ 
tion  of  the  so-called  class  matrix^]  and  inverse  halfton¬ 
ing  algorithms  for  dot  diffused  images  is  proposed  in 
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[6].  A  mathematical  description  of  dot  diffusion  algo¬ 
rithm  is  also  given  in  [7].  In  this  paper,  the  description 
of  dot  diffusion  is  reviewed  in  Sec.  2.  The  optimiza¬ 
tion  without  the  enhancement  step  is  discussed  in  Sec. 
3.  Then  the  problems  with  the  enhancement  step  will 
be  pointed  out  and  the  16  x  16  class  matrix  is  opti¬ 
mized  for  dot  diffusion  without  enhancement  in  Sec. 
3.1. 


2.  Dot  Diffusion 

The  dot  diffusion  method  for  halftoning  has  only  one 
design  parameter,  called  class  matrix  C.  It  determines 
the  order  in  which  the  pixels  are  halftoned.  Thus,  the 
pixel  positions  (ni,  112)  of  an  image  are  divided  into  IJ 
classes  according  to  (ni  mod  /,  n2  mod  J)  where  I  and 
J  are  constant  integers.  Let  x(ni,n2)  be  the  contone 
image  with  pixel  values  in  the  normalized  range  (0, 1|. 
Starting  from  class  fc  =  1,  we  process  the  pixels  for 
increasing  values  of  k.  For  a  fixed  A;,  we  take  all  pixel 
locations  (711,712)  belonging  to  class  k  and  define  the 
halftone  pixels  to  be 

h{nun2)^^Q  ifx(ni,n2)<0.5 

We  also  define  the  error  e(ni,  712)  =  x(7ii,7i2)-/i{7ii,7i2). 
We  then  look  at  the  eight  neighbors  of  (tii,  712)  and  re¬ 
place  the  contone  pixel  with  an  adjusted  version  for 
those  neighbors  which  have  a  higher  class  number  (i.e., 
those  neighbdrs  that  have  not  been  halftoned  yet). 
To  be  specific,  neighbors  with  higher  class  numbers 
are  replaced  with  a:(£,y)-h2e(7ii,7i2)/tiy  for  orthogonal 
neighbors  and  x{iyj)  +  e{ni^n2)/w  for  diagonal  neigh¬ 
bors  where  w  is  such  that  the  sum  of  errors  added  to 
all  the  neighbors  is  exactly  6(711,712).  The  extra  factor 
of  two  for  orthogonal  neighbors  (i.e.,  vertically  and 
horizontally  adjacent  neighbors)  is  because  vertically 
or  horizontally  oriented  error  patterns  are  more  per¬ 
ceptible  than  diagonal  patterns.  The  contone  pixels 
x(ni,7i2)  which  have  the  next  class  number  fc  -h  1  are 
then  similarly  processed.  The  pixel  values  a:(7ii,7i2) 
are  of  course  not  the  original  contone  values  but  the 
adjusted  values  according  to  earlier  diffusion  steps  (2). 
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When  the  algorithm  terminates,  the  signal  h(nx,n2) 
is  the  desired  halftone. 

Usually  an  image  is  enhanced  [4]  before  dot  diffu¬ 
sion  is  applied.  For  this  the  continuous  image  pixels 
C{i,j)  are  replaced  by  C'(i,  j)  =  where 

C(i,j)  =  - .  Here  the  parameter  o 

determines  the  degree  of  enhancement.  If  a  =  0|  there 
is  no  enhancement,  and  the  enhancement  increases  as 
a  increases.  If  a  =  0.9  then  the  enhancement  filter 
can  be  further  simplified[4]. 

3.  Optimization  of  the  Class  Matrix 

Knuth  introduced  the  notion  of  barons  and  near-baron 
in  the  selection  of  his  class  matrix.  A  baron  has  only 
low-class  neighbors,  and  a  near-baron  has  one  high 
class  neighbor.  The  quantization  error  at  a  baron 
cannot  be  distributed  to  neighbors,  and  the  error  at 
a  near-baron  can  be  distributed  to  only  one  neigh¬ 
bor.  Knuth’s  idea  was  that  the  number  of  barons  and 
near-barons  should  therefore  be  minimized.  He  ex¬ 
hibited  a  class  matrix  with  two  barons  and  two  near¬ 
barons.  The  quality  of  the  resulting  halftones  still 
exhibits  periodic  patterns  similar  to  ordered  dither 
methods  (See  Fig.  1).  Knuth  has  also  produced  a 
class  matrix  with  one  baron  and  near-baron,  but  un- 
fortimately  these  were  vertically  lined  up  to  produce 
objectionable  visual  artifacts.  In  our  experience,  the 
baron/near-baron  criterion  does  not  appear  to  be  the 
right  choice  for  optimization. 


Figure  I:  Dot  Diffusion  with  Knuth’s  class  matrix  (8x8 
class  matrix). 

We  used  a  cMerent  cost  fimction  for  differentiating  the 


good  halftones  from  other  halftones  [5].  In  this  cost 
function,  the  Human  Visual  System  (HVS)  is  taken 
into  account.  The  images  are  passed  through  the  HVS 
function  imitating  the  Human  Visual  System.  Since 
our  model  is  linear,  we  apply  the  HVS  function  to 
the  difference  image  between  the  original  and  halftone 
image.  Energy  of  the  resulting  image  is  defined  to  be 
the  cost  of  the  halftone  image. 

The  HVS  function  has  been  derived  in  |8]  and  in  [2]  ex¬ 
perimentally.  In  the  frequency  domain  the  HVS  func¬ 
tion  is  approximated  well  by: 


where  a  =  131.6,  b  =  0.3188,  c  =  0.525,  d  =  3.91. 
We  used  L  =  0.091  in  our  experiments  where  L  is  the 
average  luminance.  Furthermore,  the  phase  depen¬ 
dent  fimction  is  defined  as  s(^)  =  ^^^cos(4<f>)  + 
where  w  =  0.7  and  (f>  =  atan(^).  With  h{x,y) 
denoting  the  inverse  Fourier  transform  of  H{u,  v),  the 
discretized  version  /i[m,Ti]  =  /i(Tm,rn)  is  used  in  the 
calculations.  In  Fig.  2,  the  HVS  function  is  shown  for 
T=0.2. 


Figure  2:  HVS  function  H(u,v)  for 
T=0.2.  The  axes  are  ^  and 


Figure  3:  Floyd-Steinberg  error  diffusion. 


In  the  optimization  we  are  looking  for  a  class  matrix 
which  minimizes  the  cost  function.  Notice  that  the 
optimization  is  equivalent  to  finding  a  combination  of 
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numbers  from  1  to  IJ  such  that  the  related  cost  is 
minimized.  Since  the  cost  function  depends  nonlin- 
early  on  its  parameters,  we  will  use  an  optimization 
procedure  to  get  the  desired  class  matrix.  The  choice 
of  the  class  matrix  that  minimizes  this  cost  function 
was  performed  using  the  pairwise  exchange  algorithm 
[1]  described  below: 

1)  Randomly  order  the  numbers  in  the  class  matrix. 

2)  List  all  possible  exchanges  of  class  numbers. 

3)  If  an  exchange  does  not  reduce  cost,  restore  the  pmr 
to  original  positions  and  proceed  to  the  next  pair. 

4)  If  an  exchange  does  reduce  cost,  keep  it  and  restart 
the  eniuneration  from  the  beginning. 

5)  Stop  searching  if  no  further  exchanges  reduce  cost. 

6)  Repeat  the  above  steps  a  fixed  number  of  times  and 
keep  the  best  class  matrix.  ^ 


Figure  4:  Dot  Diffusion  with  the  optimized  class  matrix 
and  with  enhancement  (8x8  class  matrix) 

Example:  The  512  x  512  continuous  tone  peppers 
image  was  halftoned  by  using  Knuth^s  class  matrix 
(Fig.  1),  and  by  the  optimized  8x8  class  matrix  (Fig. 
4).  It  is  clear  that  the  dot  diffusion  method  with  the 
optimized  8x8  class  matrix  is  visually  superior  to  dot 
diffusion  method  with  Knuth’s  class  matrix.  In  fact, 
dot  diffusion  with  the  optimized  8x8  class  matrix 
offers  a  quality  comparable  to  Floyd-Steinberg  error 
diffusion  method  (Fig.  3).  Error  diffused  images  suffer 
from  worm-like  patterns  which  are  not  in  the  original 
image,  whereas  dot  diffused  halftones  do  not  contain 
these  artifacts.  Notice  that  the  artificial  periodic  pat- 

^Note  that  purwise  exchange  algorithm  yields  a  local  mini¬ 
mum  of  the  cost  function.  We  start  the  pairwise  exchange  with 
random  class  matrices  and  take  the  class  matrix  having  the  least 
local  minimum  in  order  to  get  closer  to  the  global  Tnintmiinn. 
Global  minimum  is  not  guaranteed.  _ _  _ _ 


terns  in  Fig.  1  are  absent  in  Fig.  3  and  in  the  dot 
diffusion  with  the  optimized  8x8  class  matrix  (Fig.4). 

3.1.  Dot  Diffusion  without  Enhancement 

If  we  compare  the  halftone  images  obtained  with  en¬ 
hancement  (Fig.  4)  and  without  enhancement  (Fig. 
5)  ,  we  can  conclude  that  the  enhancement  step  re¬ 
duces  halftoning  noise,  but  it  might  be  objectionable 
in  some  applications  because  of  its  very  visible  sharp¬ 
ening  effect  (e.g.,  see  Fig.  4). 


Figure  5:  Dot  Diffusion  with  the  optimized  class  matrix 
and  no  enhancement  (8x8  class  matrix). 

It  turns  out  that  we  can  get  good  halftones  without 
use  of  the  enhancement  step  provided  we  make  the 
class  matrix  larger  than  the  standard  8x8  size.  The 
price  paid  for  the  larger  class  matrix  is  that  the  par¬ 
allelism  of  the  algorithm  is  compromised.  We  found 
that  if  a  16  X  16  matrix  is  used,  the  halftone  images 
resulting  from  the  optimization  of  this  matrix  are  very 
good  even  without  the  enhancement  step.  (For  com¬ 
parison  we  note  here  that  whenever  enhancement  is 
used,  the  class  matrix  can  be  as  small  as  5  x  5  without 
creating  noticeable  periodicity  patterns.)  Such  opti¬ 
mization  was  carried  out  using  a  gray  scale  ramp  ais 
the  training  image.  The  HVS  frmction  was  used  in  the 
optimization,  and  the  associated  cost  was  optimized 
using  the  pairwise  exchange  algorithm.  The  16  x  16 
optimized  cl^tss  matrix  is  shown  in  Table  1. 

The  peppers  image  halftoned  with  the  resulting 
.  .class  matrix  is  shown  in  JFig.  6.  There  are  no  periodic- 
artifacts  in  this  result.  While  the  overall  visible  noise 
level  appears  to  be  higher  than  for  error  diffusion,  the 
problematic  halftone  patterns  of  error  diffusion  in  the 
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mid  gray  level  are  eliminated  here.  (Examine  the  body 
of  the  middle  pepper  in  Fig.  3).  By  comparing  Fig. 
I  and  6  we  see  that  16  x  16  dot  diffusion  without 
enhancement  is  also  superior  to  8  x  8  enhanced  dot 
diffusion  using  Knuth’s  matrix  because  there  are  no 
noticeable  periodic  patterns  any  more,  and  there  are 
no  enhancement  artifacts. 


Figure  6:  Dot  Diffusion  with  HVS  optimized  16x16  class 
matrix  and  no  enhancement. 

4.  Conclusion 


Dot  Diffusion  offers  more  parallelism  than  error  dif¬ 
fusion  and  the  method  has  been  optimized  in  order 
to  get  rid  of  the  periodic  artifacts.  The  enhancement 
step  prior  to  dot  diffusion  was  preserved  in  previous 
optimizations.  Since  the  enhancement  can  be  objec¬ 
tionable  in  some  cases,  the  method  has  been  improved 
by  optimizing  a  bigger  class  matrix. 
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